Small scale anisotropy in turbulent shearless mixing 
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The generation of small-scale anisotropy in turbulent shearless mixing is numerically investigated. 
Data from direct numerical simulations at Taylor Reynolds' numbers between 45 and 150 show that 
there is not only a significant departure of the longitudinal velocity derivative moments from the 
values found in homogeneous and isotropic turbulence, but that the variation of skewness has an 
opposite sign for the components across the mixing layer and parallel to it. The anisotropy induced 
by the presence of a kinetic energy gradient has a very different pattern from the one generated by 
an homogeneous shear. The transversal derivative moments in the mixing are in fact found to be 
very small, which highlights that smallness of the transversal moments is not a sufficient condition 
for isotropy. 

PACS numbers: 47.27+, 47.51+ 
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Laboratory and numerical results are continuously 
being generated on the small-scale features of turbu- 
lence dynamics (see e.g. (J-H)]). Turbulent flows con- 
tain a wide range of scales and each range is char- 
acterized by its own physics. For example, energy 
dissipation takes place at small-scales. However, the 
process is linked to the large scales of the system and 
the existence of a long-range interaction should not 
be excluded. 

A manifestation of the non-universal behavior 
of small-scales is closely related to small-scale 
anisotropy, which can be represented in terms of ve- 
locity derivative statistics. In this paper, we have 
analyzed the velocity derivative statistics of shear- 
less turbulent mixing in temporal decay, a very mild 
instance of inhomogeneity where near-isotropy could 
reasonably be expected. In fact, the turbulent shear- 
less mixing layer is possibly the simplest instance of 
an inhomogeneous turbulent flow, because it is gen- 
erated by the interaction of two isotropic turbulences 
in the absence of a mean shear flow. Therefore, in 
this flow, there is no production of turbulent kinetic 
energy and no mean convection. In general, this mix- 
ing shows that the behavior of statistical turbulent 
quantities are influenced by the presence of a kinetic 
energy gradient. 

The shearless mixing is a flow where a signifi- 
cant level of anisotropy is observed at large and 
small scales. The anisotropy persists at the moderate 
Reynolds numbers that have been reached in the nu- 
merical simulations ||-[l2j and at the moderate/high 
numbers that have been reached in the laboratory 
[l2l [l3j . These studies show that the one-point ve- 
locity statistics exhibit high intcrmittcncy in the ve- 
locity component along the mixing, as indicated by 
the large maxima of skewness and kurtosis, and only 
a mild anisotropy of the second order moments. The 
level of intermittency is a function of both the energy 



and integral scale gradients but a kinetic energy vari- 
ation alone is sufficient for the onset of intermittency 
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It should be noted, that in this flow the turbulence 
structure is different from the homogeneous sheared 
turbulence [THl, H| , but also from the turbulence near 
the fluctuating interfaces at the outer edges of turbu- 
lent shear layers, with and without free-stream tur- 
bulence HI - lij]. Here, the turbulent energy gradient 
is imposed and cannot be intensified by the contin- 
uous fluctuating interfaces produced by the instabil- 
ity of the mean shear. However, inside the shear-free 
mixing a front of high intermittency is produced and 
is displaced towards the low energy side of the flow 
@) El EH • The entrainment process is active and is 
carried out at the level of both the large and the small 
fluctuations. 

Only one inhomogeneous direction is present in this 
flow configuration. The correlations are axisymmetric 
and only the longitudinal derivative moments are sig- 
nificantly different from zero. This kind of anisotropy 
is an intermediate situation between isotropic turbu- 
lence (where only longitudinal moments are present 
and equal, see, for example, [H [T3, EH ) > an d homo- 
geneous sheared turbulence, where both longitudinal 
and transversal moments are generated, see |2|, In 
the last case the literature shows that, as the Reynolds 
number is increased, the longitudinal derivative mo- 
ment in the direction of the mean flow increases while 
the transversal odd moments (of order 3,5 and 7) are 
small in comparison to the even moments. It should 
be noticed that the information on the other two lon- 
gitudinal derivative moments (in a Cartesian reference 
frame, those in the plane orthogonal to the mean flow) 
is not available at the moment. 

We have analyzed data from numerical simulations 
in which Navier-Stokes equations have been solved in 
a parallelepiped domain with a fully dealiased (3/2 
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FIG. 1: Turbulent shearless mixing layer: (a) Scheme of 
the computational domain and boundary conditions, (b) 
Time evolution of the mixing layer thickness, varying the 
initial energy ratio Ei/E?. The mixing layer thickness 
A(f) is conventionally defined as the distance between the 
points with normalized energy values [E — E\)/(Ei — E2) 
equal to 0.25 and 0.75. The color band is the power law 
fitting with exponent 0.46. Exponent 0.33, which is indi- 
cated in the figure, is the value measured after the initial 
transient of the simulations has elapsed. 



- rule) Fouricr-Galerkin pseudospectral spatial dis- 
cretization and a fourth-order explicit Rungc-Kutta 
time integration fllij . The initial conditions have been 
obtained from a linear matching of the two homoge- 
neous and isotropic fields over a narrow region - as 
large as the flow integral scale - by means of a weight- 
ing function see figure 1 a,b. As a consequence, 
the initial statistics are closer to those present in a 
homogeneous flow than to those that would emerge 
across the mixing as the two turbulent flows interact. 
A sketch of the computational domain is shown in 
figure 1 a, x is the coordinate in the inhomogencous 
direction and y\, 1/2 arc in the homogeneous direc- 
tions normal to x. The domain has an aspect ratio 
of 2 (dimcnsionless size 2-7T 2 x Air) and the resolution 
is 128 2 x 256 at Re x = 45, 192 2 x 384 at Re x = 71, 
600 2 x 1200 at Re x = 150. All the statistical proper- 
ties of the flow are computed as spatial averages over 



planes at constant x and are only a function of (x, t). 

The two neighbouring turbulent fields are isotropic 
and each field is statistically defined by the turbulent 
kinetic energy and the integral scale (or the dissipa- 
tion rate). Let us call E\ and E 2 the turbulent ki- 
netic energy per unit mass in the two isotropic regions 
(Ex > E2) and £\, £2 their integral scales. Three pa- 
rameters characterize the mixing: the Reynolds num- 
ber, the E1/E2 ratio and the £1/^2 ratio. It has been 
found that the intermittency level and the depth of 
penetration by the eddies from the high-energy re- 
gion increase when the energy and lengthscale gradi- 
ents are concordant and decrease when they are op- 
posite. Therefore, the most efficient mixing process 
takes place when the spectra of the two mixed fields 
differ in the lowest wavenumbers, see (ll| . 

In order to have only one source of anisotropy, the 
turbulent kinetic energy gradient, two sets of simula- 
tions have been considered, both with a uniform inte- 
gral scale (€1/^2=1)- The first set of simulations has 
an energy ratio fixed at 6.6 and Reynolds numbers, 
based on the Taylor microscale, of 45, 70 and 150. In 
the second set, the energy ratio ranges from 6.6 to 10 4 
while the Reynolds number is kept equal to Re\ = 45. 

The discussion focuses on the normalized third and 
fourth order one-point moments of the longitudinal 
velocity derivative, that is, on the skewness and kur- 
tosis. These are defined as 



San/a* = {du/dxY/{{du/dxYfi 2 , 
WIW/(WJW) 3/2 , 
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K dv/dy = {dv/dyf/iidv/dyf) 
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where y is any direction normal to x (y = yi, 2/2)- 
The overbar denotes the statistical average, which has 
been approximated by a spatial average on the planes 
parallel to and inside the mixing layer at a constant 
x, see fig.l a. The turbulence is homogeneous in these 
planes. The velocity fluctuation u in equations ([TJ [2]) 
is the component of the velocity vector that is respon- 
sible for the energy transport across the mixing. 

The spatial distribution of the longitudinal deriva- 
tive skewness and kurtosis across the mixing layer is 
shown in figure [2] at several times, where two mix- 
ings with the same energy ratio of 6.6, but a dif- 
ferent Reynolds number, are compared. The spa- 
tial coordinate x has been rescaled with the mixing 
layer thickness A(i), conventionally denned as the 
distance between the points with normalized energy 
(E-E 2 )/(E 1 -E 2 ) equal to 1/4 and 3/4 jEEllil. In- 
side the mixing layer, negative values of 77 = x/A cor- 
respond to the highest homogeneous energy flow, pos- 
itive values to the lowest homogeneous energy flow. 
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r) = x/A(t) r) = x/A(t) 



FIG. 2: Longitudinal derivative skewness and kurtosis distributions: the dash-dot lines represent the statistics of du/dx, 
the continuous lines represent the statistics of dv/dy. The high energy region has a Taylor microscale Reynolds number 
equal to 45 in parts (a) and (c) and 150 in parts (b) and (d). The thick horizontal lines represent the values of the 
longitudinal skewness and kurtosis in the isotropic flow regions. 



It can observed that all these statistics depart from 
the isotropic turbulence value, -0.5, shown in the fig- 
ure by the horizontal line. However, the longitudinal 
derivatives exhibit different behavior in different di- 
rections. The main feature of these distributions is 
that the skewness departs from the opposite sign of 
the isotropic value: it is negative in the direction nor- 
mal to the mixing layer (we can observe values as high 
as -1.1), and is positive in the direction parallel to the 
layer (values as high as -0.05). 

The longitudinal kurtosis shows a maximum in the 
same zone where the skewness departures are ob- 
served, which is always in the low energy side of the 
mixing. The values of those peaks, for both the skew- 
ness and the kurtosis, are almost constant after the 
initial transient of the simulation, which lasts about 

1 11 

3-4 initial eddy turnover times r = l\jE-^ when 
Re\ = 45, but only about one eddy turnover time 



when Re\ = 150; here the distributions also show 
a fairly good collapse when the spatial coordinate x 
is rescaled with the mixing layer thickness A(t). A 
qualitative scheme of the behavior of the longitudi- 
nal skewness is shown in figure (3Ja). With respect to 
the isotropic situation, the mixing process produces a 
further compression of the filaments lying across the 
mixing layer and a reduction in the filaments compres- 
sion in the normal directions. All the mixings follow 
this common pattern, but the relative values of the de- 
viations from isotropy ASg u /g x and ASq v /q v depend 
on the flow parameters. When the Reynolds number 
increases, ASq v /q v and ASg v /g y / ASq u /q x decrease 
(figE^b)). An opposite behavior is seen when the en- 
ergy ratio increases, see again in^b). For large values 
of E1/E2, the mixing approaches a situation where a 
turbulent flow diffuses in a region of relatively still 
fluid, and the main effect on small scales in this limit 
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FIG. 4: Normalized probability density function of the 
longitudinal derivatives; data at t/r = 3.5 in the centre 
of the mixing layer from the simulation at Re\ = 150, 
E1/E2 = 6.6; £ = dui/dxi with i = x, yi and yi and 
erg is its root mean square. The probability density func- 
tion is computed by using 600 2 x 24 grid points in the 
mixing layer and 600 2 x 120 grid points in the high en- 
ergy homogeneous region. In the homogeneous region the 
probability density function compare well with the data 
by Ishihara et al. (figure 5 of [l7| and figure 4 of 20]), 
interpolated at the same Reynolds number, which use a 
larger statistical sample (512 3 ). 



FIG. 3: (a) Scheme of the general behaviour of the longitu- 
dinal derivative skewness in the shearless mixing layer, (b) 
Anisotropy of the longitudinal derivative statistics varia- 
tions inside the mixing layer as a function of the energy ra- 
tio at Re\ — 45 (circles) and as a function of the Reynolds 
number at E1/E2 = 6.6 (squares). All the quantities have 
been computed in the centre of the mixing layer, A5* is 
the modulus of the difference between the values of the 
velocity derivative skewness in the center of the layer and 
in isotropic condition. 



seems to be an additional negative stretching in the 
direction of the energy flow. Figure |4] shows the prob- 
ability density function of the longitudinal derivatives 
in the mixing, and highlights the longer negative tail 
of the probability density of the derivative in the ki- 
netic energy flow direction. 

A comparison of the longitudinal moments inside 
the mixing layer with the ones measured in the two 
homogeneous flows (HIT and homogeneous shear tur- 
bulence) is shown in figure [5] The data for the 
isotropic turbulence in this figure are taken from the 
homogeneous regions of the present simulations and 
from the reviews by Sreenivasan and Antonia 
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FIG. 5: Comparison of the derivative skewness modulus in 
the shearless mixingand in the homogeneous shear flow 
(data from 0] and [3]). The data for the homogeneous 
and isotropic turbulence are taken from the homogeneous 
regions of present simulations and from [lj . Data from ref . 
have been read from published graphs. 
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while the data for the homogeneous shear flows are 
taken from Warhaft and Shen 0, E| ■ These last data 
are laboratory data which, due to the high level of 
technical difficulties which characterize this kind of 
measurements, only give the longitudinal derivative 
of the velocity component in the streamwise direc- 
tion. The difference with the homogeneous shear 
flow, which generates large transversal moments, but 
has less influence on the longitudinal ones (with re- 
spect to the isotropic values), is immediately appar- 
ent. Shear and shear less flows have thus different 
kinds of anisotropy in the small scales: a strong dif- 
ferentiation of longitudinal derivative moments for 
shearless flows and high values of transversal deriva- 
tive moments for shear flows. 

In conclusion, the simulations we have carried out 
show that a significant small scale anisotropy and in- 
termittency is generated in a decaying shearless tur- 
bulent mixing. This intermittency is characterized by 
a large departure of the longitudinal derivative mo- 
ments, which are different in the directions across and 
parallel to the layer, from the typical values of the 
isotropic condition, even in this flow where there is 
no energy production (due to the lack of mean flow 
gradients). The deviations from the isotropic values 
are large and follow a common trend: the longitudinal 
derivatives in the energy gradient direction are more 
intermittent, while the intermittency is milder in the 
orthogonal directions. We also observe that a small 
intermittency on transversal velocity derivatives does 
not necessarily mean a tendency towards isotropy. 
The structure of the anisotropy is such that the skew- 
ness departure from isotropy reduces the compression 
on fluid filaments parallel to the mixing layer and 
enhances that of the filaments orthogonal to it. In 
the shear-free mixing, the small-scale turbulence has 
a different structure than in the homogeneous shear 
flow case [3-11, H| . The symmetry is different and the 
anisotropy here is due to the inhomogeneity. We think 
that this is the principal agent for the difference in 
the anisotropy structure of the small-scale in these 
two flows. The reduction of the skewness negativity 
in directions parallel to the mixing (relative elonga- 
tion with respect to the isotropic situation) and the 
enhancement of the negativity across the mixing (rel- 
ative compression with respect to isotropy) should 
also be linked to the incompressibility of the flow. 



These effects also persist at moderate Reynolds num- 
bers, therefore the asymptotic approach to the local 
isotropy, if present, is very slow. 
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